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Graphic Realizations of Joint-Degree Matrices* 
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Abstract 

In this paper we introduce extensions and modifications of the classical degree se¬ 
quence graphic realization problem studied by Erdos-Gallai and Havel-Hakimi, as well 
as of the corresponding connected graphic realization version. We define the joint-degree 
matrix graphic (resp. connected graphic) realization problem, where in addition to the 
degree sequence, the exact number of desired edges between vertices of different degree 
classes is also specified. We give necessary and sufficient conditions, and polynomial time 
decision and construction algorithms for the graphic and connected graphic realization 
problems. These problems arise naturally in the current topic of graph modeling for 
complex networks. From the technical point of view, the joint-degree matrix realization 
algorithm is straightforward. However, the connected joint-degree matrix realization 
algorithm involves a novel recursive search of suitable local graph modifications. Also, 
we outline directions for further work of both theoretical and practical interest. In par¬ 
ticular, we give a Markov chain which converges to the uniform distribution over all 
realizations. We show that the underline state space is connected, and we leave the 
question of the mixing rate open. 


1 Introduction 

Let di > d 2 dn he a. sequence of integers. The classical graphic realization problem 

asks if there a simple graph on n vertices whose degrees are exactly di > d 2 > ■ ■ ■ > dn- Erdos 
and Gallai showed that the natural necessary conditions for graphic realizability, namely that 
each subset of the highest k degree vertices can absorb their degrees within their subset and 
the degrees of the remaining vertices: Y^i=i — ^(^~l) + Dr=A;+i dj}, are also sufficient 

di- The well known Havel-Hakimi algorithm [HI 120] achieves a realization in an efficient 
greedy way. It repeatedly sorts the vertices according to residual unsatished degree, picks 
any vertex of residual degree A, and connects it to the A vertices of highest residual degree. 
The process is repeated until all the degrees are satished. If one further wants to construct a 
connected graphic realization (a requirement which is clearly important in networking), Erdos 
and Gallai showed that the obvious necessary condition — 2(n —1) (i.e., there is a 
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spanning tree) is also snfficient. In particnlar, it is easy to see that a non-connected realization 
can be transformed to a connected realization by a seqnence of flips, each flip breaking a cycle 
inside a connected component, and redncing the nnmber of connected components by one. A 
“flip” picks two edges xy and uv snch that xu and yv are not edges, removes xy and uv from 
the graph, and adds xu and yv to the graph. It is clear that flips do not change the degrees 
of the graph. 

Now, let V = [n] be a set of vertices. Let V = {Vi, V 2 ,..., 14} be a partition of V denoting 
snbsets of vertices with the same degree and let d : V —)■ N be a fnnction denoting the degree 
of vertices in class V}. Let D = (dij) he a k x k matrix denoting the nnmber of edges between 
Vi and Vj] if i = j it is the nnmber of edges entirely within V). The joint-degree matrix graphic 
realization problem is, given {Y,d,D), decide whether there is a simple graph G on V, snch 
that, Wi each vertex in Vi has degree d{Vi), Wi 7 ^ j there are exactly dij edges between V) and 
Vj, and, Wi, there are exactly da edges entirely inside V). The joint-degree matrix connected 
graphic realization problem is to decide whether a connected graphic realization for (V, d, D) 
exists. Fnrthermore, we want to either constrnct snch a realization, or ontpnt a certificate that 
no snch graphic connected realization exists. In this paper we give necessary and snfficient 
conditions, and polynomial (in n) time algorithms for the decision and constrnction of the 
joint-degree graphic and connected graphic realization problems. 

The practical significance of joint-degree matrix realization problems arises in graph mod¬ 
els for several classes of complex networks. For instance, in networking, models for Internet 
topologies are constantly nsed to simnlate network protocols and predict network evolntion. 
Commonly nsed topology generators, snch as GT-ITM |3H [ 6 ] which generates random nearly 
regnlar graphs, and power-law based generators uniiaiiiiETiEsmn], generate random graphs. 
However, the properties of topologies constrncted nsing pnrely random graph models were 
challenged, most notably in [Sj El]. Using qnalitative argnments and striking images, they 
argned that power law random graphs constrnct a dense core of nodes with very high degrees, 
while nodes of smaller degree are mostly attached to the periphery of the network. On the 
other hand, highly optimized Internet topologies place low degree bnt very high bandwidth 
renters at the center of the network, while high degree nodes are mostly placed in the pe¬ 
riphery to split the signal manyways toward the end nsers. To qnantify their argnment, [21| 
used a random graph Gr(U,Ur) under several power law models [U EH [T] [301 EHi a 
real network topology G(y,E). They found that deg(M)deg(n) is much larger than 

deg(M)deg(n). Independently, [28l ES] made the same observation for several other 
technological and biological networks. 

Going one step further, E5] argued that, a determining metric for a graph of given 
degrees to resemble a real network topology, is the specific number of links between vertices 
in different degree classes. Using heuristics that presumably approximate the target number 
of edges between degree classes, [25| constructed graphs strikingly similar to real network 
topologies. The joint-degree matrix graphic and connected graphic realization problems stud¬ 
ied in Sections and formalize the approach of |2Sl [2S] • In general, one would want to 
construct a uniformly random realization of {Y,d,D). We state this as an open problem in 
Section 4. However, for practical purposes, the heuristics of [26l EH] achieved very satisfactory 
results using randomness in a configuration model adjusted to the (V, d, D) problem. On the 
other hand, there are no theoretical results concerning the properties of this model. In the 
same sense, our joint degree matrix realization algorithm in Section 2, also allows substantial 
randomness in the choice of the edge to be added in each of its greedy steps. 
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In Section we address the joint-degree matrix graphic realizability problem. We show 
that the natural necessary conditions are also sufficient, and can be checked efficiently. We 
also obtain a second polynomial time construction algorithm, which is a key for the algorithm 
in Section This Balanced Degree Algorithm constructs the graph by increasing the number 
of edges, without increasing the number of connected components. In Section we address 
the joint-degree matrix connected graphic realizability problem. By sharp contrast to the 
degree sequence connected realization, here the necessary and sufficient conditions are fairly 
complex, and of exponential size. However, using a recursive algorithm that searches for suit¬ 
able local graph modihcations to construct a connected graph, we manage to either construct 
such a graph in polynomial time, or identify at least one necessary condition that fails to 
hold. In Section]^ we discuss structural differences between degree sequence and joint-degree 
matrix problems. In particular, the former are known to be related to matchings, while no 
corresponding fact is known for the latter. Finally, in Section]^ we propose a natural Markov 
chain for sampling from (V, d, D), and we show that it is ergodic. 

Recent related work. Independently, |32l HOl HE] give polynomial time algorithms for 
constructing a graph in {Y,d,D). Moreover, in |32] an alternative proof is proposed for the 
fact that the Markov chain we dehne in Section is ergodic. This proof however was flawed, 
as noted in uni, where an alternative proof is given as well. With respect to the mixing time 
of this Markov chain, [32] performed experiments based on the autocorrelation of each edge; 
these experiments suggest that the Markov chain mixes quickly. In a more recent work, [T5] 
shows fast mixing for a related Markov chain over the subset of (V, d, D) that contains the 
balanced realizations, i.e. realizations where for each i,j the edges connecting Vj to Vj are 
as uniformly distributed on Vi as possible. (Notice that this is not what we call a balanced 
graph, e.g. in Lemmaor in Section]^) 


2 Joint-Degree Matrix Graphic Realization 

Let V = [n] be a set of vertices and V= {Hi, V 2 ,..., 14} be a partition of V denoting subsets 
of vertices with the same degree. Let d : V —>■ N be a function denoting the degrees of vertices 
in class H and D = (dij) he a k x k matrix denoting the number of edges between H and Vj] if 
i=j it is the number of edges entirely within H- The joint-degree matrix graphic realization 
problem is, given (V, d, D), decide whether there is a simple graph G on V, such that, Vi each 
vertex in V} has degree d(H), Vi 7 ^ j there are exactly dij edges between H and Vj, and Vi, 
there are exactly da edges entirely inside H- We use the notation (V, d, D) to also denote the 
set of all such graphs. 

We will prove that the following natural necessary conditions for the instance (V, d, D) to 
have a graphic realization are also sufficient: 

(i) Degree feasibility: 2dii + J2j(zik]j^idij = |H| ■ d{V), for 1 < i < fc. 

(a) Matrix feasibility: The matrix D is symmetric with nonnegative integral entries, and 
dij < I H| ■ for 1 < i < j < /c, while da < |H| ■ (|Ti| —1)/2, for 1 < i < /c. 

There is a straightforward algorithm for constructing a graph G G {Y,d,D). First, the 
algorithm constructs a graph G' that has the “right” number of edges between any Vi, Vj 
(or within a Vj)- Then, the degrees within each H are taken care of, resulting in a graph 
G e {Y,d,D). 

The algorithm proceeds as follows: 
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Start with an empty graph G on V 
For each i 

choose arbitrarily da edges between vertices of and add them to G 
For each j 

choose arbitrarily dij edges between vertices of Vi and and add them to G 
For each i 

While not all degrees in Vi are equal 

Choose u,v EVi such that deg(M) < diV) and deg(n) > diyi) 

Find min{|deg(M) — d(ld)|, |deg(n) — diVi)]} neighbors of v that are not neighbors of u 
Disconnect them from v and connect them to u 
Output G 

To see that the algorithm works, first notice that if dij < |Di| ■ |V^|, for 1 < i < j < k, and 
da < |fd| ■ (|fd| —1)/2, for 1 < i < /c the edge-adding part of the algorithm works. This results 
in a graph G that satishes the dij requirements, but not necessarily the degree requirements. 

Now, assume that there exist some i such that not all the degrees in Vi are equal to diVi). 
If ‘2dii+J2j(zikij^i dij = |Id| ■d(Id), this means that there exist u,v E Vi such that deg(M) < d{V) 
and deg(n) > diVi). Also, there are |deg(M) — deg(n)| neighbors of v that are non-neighbors 
of u, and |deg(M) — deg(n)| > min{|deg(M) — |deg(n) — d{yi)\}. Also, notice that each 

iteration in the “while” loop reduces the number of “wrong degrees” by at least one, without 
affecting the dij requirements. That is, in at most n iterations G E {Y,d,D). 

Although the Joint-Degree Matrix Graphic Realization problem has a straightforward 
solution, this is not the case if we also ask for the resulting graph to be connected. Before 
we move to this problem, we present an alternative algorithm for the Joint-Degree Matrix 
Graphic Realization that we will need later. 

Balanced Degree Algorithm 

This construction algorithm grows the graph G in iterations, one edge at a time, starting from 
the empty graph Gq, keeping the edges between each Vi and Vj always < dij (resp. the edges 
within each Vi always < da), and ending with a realization G E (V, d, D). The key idea of the 
algorithm is to maintain a balanced degree invariant within each Vi. If is the graph after 
iteration i, the algorithm maintains max^gv^ degc.^(n) — min^gy. deg(j^(n) < 1, for 1 < i < A; 
(where degQ(n) is the degree of vertex v in graph G, as usual). This motivates the following 
dehnition. For the graph Gi after iteration £, for all 1 < i < /c, let A^i,G£ = {v E Vi ■. degQ(n) = 
min^gy. degG.^(M)} and let Mi^Gt = {v E V : degG.^(n) = max^gy. degG,(M)}. 

The algorithm proceeds as follows: While there is some i and j (possibly i = j) such 
that dij is not satished, the construction algorithm picks any such i and j and adds an edge 
between Vi and Vj (resp. inside V), while maintaining the balanced degree invariant and 
without affecting the extend to which the other (i„„’s are satisfied. Let us assume that we 
are at the beginning of the (£-|-l)th iteration, and i,j have been picked such that dij is not 
satished. Let G = Gi and assume W = W and Mi = Mi^Ci suitably dehned for all i. 
There are several cases to consider: 

If i ^ j consider Cases Al, A2 and A3 below, in the order that they are listed: 

Case Al: if there exist u E Ni,v E Nj such that uv ^ E{G) then add uv to E{G); 

Case A2 if there exist u E Ni,v E Mj such that uv ^ E{G) then 
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pick a v' E Nj and find a neighbor x of v such that v'x ^ E{G)] 
delete the edge vx from E{G) and add the edges uv and v'x to E{G); 

Case A2': if there exist u G Mj,n G Nj such that uv ^ E{G) then symmetric to Case A2; 
Case A3: hnd u G Mi,v G Mj such that uv ^ E{G); 

pick u' G Ni and hnd a neighbor x of u such that u'x ^ E{G)] 
pick v' G Nj and hnd a neighbor y of v such that v'y ^ E{G)] 
delete the edges ux,vy from E{G) and add the edges u'x,uv,u'y to E{G); 

If i = j consider Cases Bl, B2 and B3 below, in the order that they are listed: 

Case Bl: if there exist u,v E Ni such that uv ^ E{G) then add uv to E{G)] 

Case B2 if there exist u G Ni,v E Mj such that uv ^ E{G) then 
if |Aj| = 1 then add uv to E{G) 
elseif \Ni\ > 1 then 

pick a v' E Ni and hnd a neighbor x of v such that v'x ^ E{G)] 
delete the edge vx from E{G) and add the edges v'x and uv to E{G); 

Case B3: hnd u,v E Mi such that uv ^ E{G)] 
if |Aj| = 1 then 

pick a w E Ni and hnd a neighbor x of u such that wx ^ E{G)] 
delete the edge ux from E{G) and add the edges wx and uv to E{G)] 
elseif |Aj| > 1 then 
pick w, w' E Ni] 

hnd a neighbor x of u such that wx ^ E{G)] 
hnd a neighbor y of v such that w'y ^ E{G); 

delete the edges ux,vy from E{G) and add the edges wx,uv,w'y to E{G)] 

Theorem 1. If the degree and matrix feasibility conditions hold, then the above algorithm 
constructs a graph G E {Y,d,D). The algorithm runs in time polynomial in n. In particular, 
\/i, if Gi is the graph at the end of the ith iteration, then, the number of edges between Vi and 
Vj (resp. inside Vi) have increased by one, the number of edges between and inside all other 
degree classes have not changed, and the balanced degree invariant holds. 

Proof. Assume that we are at the beginning of the f'th iteration and the balanced degree 
invariant holds. We show that the above algorithm maintains the invariant after the next 
edge is added. 

Observe hrst that, for all i, the sets N^ and Mj are always nonempty and Mj U Aj = V^. 
In fact, either Aj = Mj = Vi, or {Mj, Aj} is a partition of Vj. In the sequel, if we refer to Mi, 
then we are assuming that Aj ^ Mi. 

Now let i and j be two indices so that Ge has less than dij edges in the subgraph induced 
hy VUVj. By matrix feasibility, there is an edge uv ^ E{Ge) with u E Vi and v E Vj. 

Consider first the case i ^ j. If the edge uv falls into Case Al, the invariant clearly holds 
for G^ + uv. In Case A2 (A2' is symmetric) since Mj ^ Nj, there is a v' ^ v E Nj. Since 
deg(n') < deg(n), there is a neighbor x of v such that v'x ^ G^. The specihed actions maintain 
the invariant. 

If Case A3 is reached, we must have uv such that u E Mi and v E Mj. Consider u' E 
Ni,v' E Nj. Since no edge in Case A2 (or A2') was available, we must have uv',u'v and 
u'v' E E{G(). Since n is a neighbor of u' but not of u, and degg(M') < deg(^(M), there must 
exist some x that is a neighbor of u but not of u' . Similarly, there must exist some y that is 
a neighbor of v but not of v' (possibly y = x). Because u, u' are both in Vj and v, v' are both 
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in Vj, we can remove xu, yv and add xu', yv'. This way, the extend to which the requirements 
of matrix D are satisfied is not affected, but the degrees of u,u',v and v' change so that we 
can add uv to G and the statement of the theorem is true. 

Next consider the case i = j. If u,v G iVj, then it is clear that Case B1 maintains the 
invariant. 

Now, assume there is no available edge with both ends in Ni, but there exist u G Mi, v ^ Ni 
such that uv ^ E{G). If |iVj| = 1, adding uv to the current graph maintains the invariant 
since after the addition Mj = {n}, Ni = Vi — {n}. 

So suppose that |iVj| > 2 , and let v' ^ u E N^. Since deg(r;) > deg(n'), there exists an edge 
vx G Gi such that v'x ^ G^. Note that x ^ u. Now the specified actions satisfy the theorem. 

The last possibility is that the only available edges have G Mj 7 ^ Ni. Notice that by 
exhausting Case B2 hrst, w E Ni implies wu,wv G E{G(). Again, we consider two cases: 
m = 1 and |iV,| > 2 . 

In the former case, pick a. w E Ni. Since n is a neighbor of w but not of u, and degciw) < 
degciu), there must exist some x that is a neighbor of u but not of w. Because u, w both are 
in Vi we can remove xu and add xw. This way the dij requirements are not affected for any 
j, but the degrees of u and w change so that we can add uv to G and keep the invariant true. 

In the latter case, where |A^i| > 2, pick a w,w' E Ni. Notice that tc is a neighbor of w' 
(or else ww' would have been added in Bl). Since n is a neighbor of w but not of u, and 
degciw) < degciu), there must exist some x that is a neighbor of u but not of w. Similarly, 
there must exist some y that is a neighbor of v but not of w' (possibly y = x). Because 
u,v,w,w' all are in I 4 , we can remove xu, yv and add xw, yw'. This way the dij requirements 
are not affected for any j, but the degrees of u,v,w and w' change so that we can add uv to 
G and the invariant holds. 

As long as there exists some dij not yet satisfied, the algorithm manages to increase 
the number of Vi-Vj edges by one without changing the number of Vii-Vji edges, for any 
7 ^ {hi}- Therefore in m = J2i<j^i,j iterations, all the edge requirements are met. 
Now, since degree feasibility holds and for any i we have 2dii + J2jjLi dij = \Vi\ ' diVi), we get 
that Vn G Vi, degciy) = diVi), as desired. □ 

Remark 1: The transformations of adding and deleting edges in all non-trivial cases of the 
algorithm resemble augmenting paths. However, in general, these transformation are not 
augmenting paths. For example, in Case A3, the sequence of edges u'x, xu, uv, vy and yv' 
includes the case where x = y. We clearly have an alternating sequence but not a path. We 
shall revisit this comment in Section [5l 

Remark 2: We claim that the construction algorithm never increases the number of connected 
components. In particular, it can be verihed that, in all cases, when an edge is removed, a 
path between its endpoints is created by the edges added in the same iteration. In particular, 
if the graph G^ at iteration i is connected, the hnal output graph G will be connected. We will 
use this fact critically in the algorithm which constructs a connected realization of (V, d, D) 
in Section m 

A generalization 

It is natural to consider the generalization of the joint degree matrix problem (V, d, D*), where 
we allow the entries of D* to be in Nq U {*}. If dij = *, there is no restriction on the number 
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of edges between the corresponding sets. We can use the main idea of the proof of Theorem 
to provide a polynomial time construction algorithm. This is proved in Theorem below. 

Notice that, if the entire matrix D* consists of =t:’s, then this is the standard degree sequence 
realizability problem. In the case where D* contains both integers and *’s, if (V, d, D*) is 
nonempty, then there exists some graph G G (V, d, D*) such that the subgraph H of G 
dehned by the integer entries of D* satishes the balanced degree invariant. This is proved in 
Lemma 1^ below. We call such a graph G a balanced {Y,d,D*) graph. 

Lemma 2. If {Y,d, D*) ^ 0, then there exists a balanced {Y,d,D*) graph. 

Proof. Let G G {Y,d,D*) and consider the subgraph H of G with vertex set V{H) = V{G) 
and edge set E{H) = {uv G E{G) : u ^ Vi, v E Vj, dij ^ *}. Assume that H does not 
satisfy the balanced degree invariant. Then consider G' on V{G) with edge set E{G') = 
E{G) \ E{H). We can hnd u and v in some Vi such that, degj:^(M) > degfj{v) + 1, and thus 
degQ/(M) < deg( 5 /(n) — 1. We can pick a neighbor x of u in H that is not a neighbor of v, and 
a neighbor ?/ of n in G' that is not a neighbor of u. We remove ux, vy in G and add vx and 
uy. We repeat the above procedure until H satishes the balanced degree invariant. Notice 
that the edge hips are such that the resulting graph G is still in (V, d,D*). □ 


Theorem 3. If {Y,d,D*) ^ 0, then we can construct a graph {Y,d,D*) in polynomial time. 

Proof. Notice that for any balanced graph G G {Y,d,D*), if we consider the subgraph H 
dehned by the * entries of D*, it is a realization of the same degree sequence, say d'l, d' 2 , ..., 
and the same edge restrictions, i.e., there are no edges between vertices of Vi and Vj if dij 7 ^ *. 
This fact, together with Lemma suggest the following algorithm to construct a graph G G 
(V, d, D*), if one exists: Let Dq be the matrix we get if we substitute all *’s with O’s. We hrst 
run the construction algorithm of this section with input Y,Dq. This will construct a graph 
H that has dij edges between Vi and Vj if dij * and 0 otherwise, and satishes the balanced 
degree invariant. Now, what we need to add to get a G G (V, d, D*) is a graph on V with 
degree sequence dehned by diV) — clegji^(M) for every u E Vi, where Vi-Vj edges are forbidden 
if dij 7 ^ *. However, this is reduced to hnding a matching in a properly dehned graph Q. 

To see this, notice that if we want to construct a graph on vertex set [n] with a given 
degree sequence di > d 2 > ... > dn and a given set E of forbidden edges, we can dehne the 
graph Q = {V,S) as follows. The vertex set is V = where Vi = ^j&[n],jf^^{vij} 

{uij/}, \/l < i < n. The vertices Vij denote a potential edge between i and j in the graphic 
realization. The vertices Wjy will enforce the required degrees di. Now the edges are S = 
dli<i<j<n, ij^Fiivij, Vji}} U]ll {{vij, Uij>}}. It is straightforward to verify that the 

degree sequence, given E, is realizable if and only if Q has a perfect matching. □ 


3 Joint-Degree Matrix Connected Realization 

We now turn to the question of constructing a connected graphic realization of an instance 
(V, d, D), or showing that such a realization does not exist. It is easy to see that this problem 
is different from its counterpart Erdds-Gallai condition (the degrees summing up to at least 
2(n —1)). In particular, there are graphically realizable instances of {Y,d,D) which include 
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many edges, but have no graphic connected realization (for example, if all the edges are 
required to be inside distinct classes Vi and Vj, i^j). It is also easy to see that arbitrary simple 
flips cannot be used to decrease the number of connected components of a non-connected 
graphic realization G G {Y,d,D). In particular, let uv and xy be edges in G, let ux and vy 
be edges not in G, and let u E Vi, v G Vj, x E Vi/ and y EVf. Then, the flip of removing uv 
and xy and adding ux and vy yields a graph in (V, d, D) if and only if = Vi/ or/and Vj = Vy. 

In what follows we give necessary and sufficient conditions for (V, d, D) to have a connected 
realization. The proof provides a polynomial (in \V\ = n) time algorithm that constructs a 
connected realization, if one exists, or produces a certificate that (V, d, D) does not have a 
connected realization. 

Roughly, the general approach to construct a connected graph in (V, d,D) is to first 
construct a tree on V that does not violate the upper bounds specified by D. This will 
be called a valid tree for (V, D). If such a tree exists, then Lemma shows how to transform 
it to a tree that does not violate the upper bounds specified by D and d and also satisfies 
the balanced degree invariant.We call such a tree a balanced tree for {Y,d,D). We may then 
continue with the greedy construction algorithm of Section which, by Remark 2 at the end 
of Section never increases the number of connected components, so that we extend the 
balanced tree to a (V, d, D) graph. 

Lemma 4. Let Tvaiid be a valid tree for (V, D). Then, we can efficiently construct a balanced 
tree for (V, d, D), Tsaianced- 

Proof. To construct Teaianced we modify the degrees within each Vi. We make use of the fact 
that, in any tree, if we pick two vertices u and v, we can move any neighbor of u to become 
a neighbor v, with the exception of the neighbor that lies on the unique u-v path of the 
tree. Let 5i be the average degree in Tvaiid of the vertices in Vi. Then, as long as 3u,v E Vi 
such that deg^T^^j.^ (m) > [5*] and degy^^j.^(u) < move neighbors of u to v, until either 
degTvaiid ('^) ~ ~ This will need at most 2|Vi| iterations. We do this 

for every i to get a balanced tree Teaianced- Notice that, while the degrees are made as equal 
as possible, the number of edges between Vi and Vj is not affected, for any i,j. Thus, Teaianced 
is a balanced tree for (V, d,D). □ 

If there is no connected realization, then we want to produce a certificate of non existence of 
a valid tree. In general, it is not clear how to construct efficiently a valid tree or a certificate of 
non existence of such a tree. Indeed, the sufficient and necessary conditions for connectivity 
listed below appear to require exponential search. However, our Valid Tree Construction 
Algorithm solves both problems in polynomial time. With this motivation in mind, we proceed 
to the technical details. 

We need some more notation. Let G be a connected realization in {Y,d,D). Let V = 
{Vi,... ,14}, where |V| = max{l, \Vi\ — da} is the minimum possible number of connected 
components of a subgraph of G induced by V- Let D = (dij) be derived from D in the natural 
way: dii = 0 and djj = min{|Vi| ■ \Vj\,dij}- 

Lemma 5. A valid tree T for (V, D) can be transformed efficiently to a valid tree T for (V, D). 

Proof. To turn T to a tree T on V, that satisfies the dj/s as upper bounds we stick a path 
of length \Vi\ — \Vi\ on an arbitrary vertex of Vi, so that we get |V| vertices. We do this for 


every i. Now notice that for the resulting tree T we have that the edges of T inside Vi are 
|hi| — \Vi\ — 1. Recall that, by dehnition, \Vi\ = max{l, \Vi\ — da}. Therefore, the edges of T 
inside Vi are at most da. Moreover, those paths inside each Vi were the only thing added to 
T. That is, the number of edges between Vi and Vj are at most dij < dij. Thus we created a 
tree T on V that does not violate the dj/s, and therefore the dj^’s. □ 

Therefore, Lemmata and reduced the problem to hnding sufficient and necessary con¬ 
ditions for the existence of a valid tree T for (V, D), an efficient construction for some T, if it 
exists, or a certihcate that T does not exist. 

We are ready to state the necessary and sufficient conditions. Let = {1^^, ... 1^^} C 

{Ri, ^ 2 , ■ ■ ■, 14}, and let A = {Mi, M 2 ,... Ma} be a partition of T. The interpretation is that 
each Mi will collapse to a single vertex. Dehne the undirected weighted graph Q = (V, T, vS) as 
follows: 

• V = (ai, 0 : 2 ,... ax} ^j^i-i,i 2 ,...,ie {uj}, that is, one vertex a* for each M* and one vertex uj for 
each Vj ^ J^. 

• If there exists I 4 G Mi and Vy G Aj such that d^y > 0, then aiaj G £ and w{ai, aj) = 1. 

• If there exists Vj ^ J^, and, for some i, there exists I 4 G M, such that d^j > 0, then UiUj G £ 
and w{ai, Uj)=mm{\Vj\,J2xVeAi 

• If there exist Vi ^ T ,Vj ^ T and dij > 0, then UiUj G £ and w{ui,Uj) = dij. 

Necessary and Sufficient Conditions: Given a connected G G (V, d, D) we can easily get a valid 
tree T for (V, D). To do so, we collapse each connected component of the subgraph Gi of G, 
induced by I^, in to a single vertex, and if necessary a few of these vertices together, so that the 
cardinality of the vertices of Gi reduces from \ Vi\ to |Vi|. Then, we delete any loops or multiple 
vertices. The resulting graph is still connected. We do this Vi and then take a spanning tree 
of the resulting connected graph. This is a valid tree T for (V, D). Then, the existence of T 
implies the following necessary condition for a connected realization to exist: for every iF and 
every M, the graph ^ = (V,T,tc) is connected and — 1*^1 + 'Ylii-Vi^T 1^*1 ~ 

On the other hand, we call a pair (J-”, M) such that the above necessary condition fails, a 
certificate that no connected realization of (V, d, D) exists. 

Next, we will show that the above stated necessary condition for a connected realization 
of {Y,d,D) to exist is also sufficient. In particular, we will prove that the Algorithm Valid 
Tree Construction below either produces a valid tree T for (V, D), or produces a certihcate 
{F, M) that no connected realization of (V, d, D) exists. 

The construction algorithm is as follows. Let V = The algorithm tries to construct 

a valid tree on V by maintaining |R|—1 edges which are valid for (V, D) (the number of edges 
between each Vi and Vj never exceeds dij), while at the same time decreasing the number of 
connected components by adding and removing edges appropriately. The main idea is that, 
if two components cannot he connected in a trivial way that maintains validity, then the Vi’s 
that intersect more than one connected components play a critical role. We constantly try to 
“free” an edge incident to such a Vi while preserving validity and not increasing the number 
of connected components. In the case that such a Vi intersects a cycle, this is an easy task. 
Otherwise, we have to remove all the Vfs that intersect more than one components, and try to 
connect two components in the resulting graph. We recursively repeat this until we connect 
something, and then it is easy to hnd a sequence of adding and removing edges that connects 
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two components in the original graph and maintains validity. If the recursion fails, we have a 
certihcate that no connected realization exists. 

Algorithm Valid Tree Construction (V, D) 
begin 

V = start by a graph Go consisting of |C| —1 valid edges over V; 

j = 0; comment: j is the depth of the recursion; 

G = Go; comment; G is an auxiliary graph; 
while Go is not connected 

begin 

Aj = {v : V lies in some cycle of Gj}; 

Cj = {Vi : Aj n Vi ^ 0}; comment: V^’s intersecting some cycle of Gj; 

Pj = {Vi : Vi intersects at least two connected components of G^}; 

Zj = {e G Gj : at least one endpoint of e is in some V E Pj}; 

Case 1: if 3 an edge e ^ Go connecting two connected components in G 
without violating upper bounds of D 
add e to G; remove any edge from any cycle of G; 
j = max{j — 1,0}; Gj = G U Pj U Zj] G = Gj; 

Case 2: elseif Gj fl Pj 7^0 

pick u,v in some Vi E Gj (1 Pj that are not connected in Go and u E Aj A Vi] 

hnd a neighbor x of u that lies on the same cycle in G ; hnd any neighbor y of a; in G; 

remove xu and yv from G and add xv and yu to G; 

j = max{j — 1,0}; Gj = G U Pj U Zj] G = Gj] 

Case 3: elseif Pj = fli 

let Gi, G 2 ,..., G^ be the connected components of Gj] 

let Ai = {Vx : VxPV{Gi)} foi l<i <^]\et P=ul^^Ai]\et A = {Ai,... ,Ai}] 
output (P, A) and terminate; comment: found a certificate of non existence] 

Case 4: else Gj+i = Gj \ Pj] G = Gj+i] j=j + 1; 

end; 

output Go; comment: found a connected realization] 

end. 

Theorem 6. Algorithm Valid Tree Construction outputs a valid tree for (V, P), if such a 
tree exists. Otherwise, it outputs a certificate (P, A) showing that no such tree exists. The 
algorithm runs in time polynomial in n. 

Proof. First, notice that we start with a graph Go on \V\ vertices and \V\ — 1 edges, that 
does not violate any upper bounds imposed by D. The algorithm also creates at most k 
graphs Gi, G 2 ,..., Gfc, such that G^+i is an induced subgraph of Gi (and thus of Go as well) 
on strictly less vertices. To be precise, V{Gi) \ V(Gj+i) contains one or more of the h}’s. To 
see this, notice that for Gj+i to be constructed this must happen in Case 4, and then Gj+i is 

the subgraph of G* induced on V{Gi) \ . But Vi in this case should be nonempty, 

or the algorithm would have terminated in Case 3. Now, for each of these graphs, say Gj, we 
associate four sets: ^ 

Vj is the set of the Vfs that intersect more than one connected component of Gj, 

Zj is the set of edges that have at least one endpoint in some Vi EVj, 
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Aj is the set of all vertices that belong to some cycle in Gj, 

Cj is the set of the Vi’s that intersect some cycle in Gj. 

Notice that if Go is not connected, then for any j such that a Gj is constructed by the 
algorithm we have Aj 7^ 0 (and thus Cj 7^ 0 as well). To see this notice that Go is either a tree, 
or contains a cycle. Moreover, whenever a Gj is created (in Case 4), its vertex set contains 
Aj_i (otherwise, the current iteration would not go further than Case 2). That is, if Gj_i 
contains a cycle, so does Gj (if created at all). The above also implies that whenever Gj is 
created, V{Gj) 7^ 0. As discussed above, the algorithm cannot go to Case 4 for k consecutive 
iterations. Therefore, within k iterations one of Cases 1,2 or 3 happens. 

Suppose that either Case 1 or Case 2 happen, i.^, two connected components G, G' of 
G = Gj become connected to each other without any dij being violated. If j = 0, the number 
of connected components of Gq is decreased. Assume not. Notice that G and G' in Gj_i must 
be subgraphs of the same connected component, otherwise they would have been connected to 
each other in an earlier iteration, when G was still Gj-i. Now, Gj_i is updated by adding Vj-i 
and Zj_i back to G and this becomes the current graph G. Just for notational convenience, 
we are going to call this graph G)_^ as opposed to the old Gj_i. Notice that Vj-i is not 
affected, i.e., 'Pj_i = Vj-i. The key observation now is that, in G)_^ a new cycle is created, 
containing one new edge added in the last iteration, as well as some v from some Vi G Vj-i 
(that was on a path connecting G and G' in Gj_i). That is, in the next iteration Case 2 will 
happen and the number of components of G)_^ will go down by one as well. So, we have that 
if Cases 1 or 2 happen, while G = Gj, then in the next j iterations Case 2 will happen. This 
results in decreasing the number of connected components of Go by one. 

Now suppose that Case 3 happens. That is, for some j, Vj = 0 and Cases 1 and 2 fail. 
Let Gi, G 2 , ..., G^ be the connected components of Gj and let A and be dehned as in the 
algorithm. The dehnition of the Mj’s makes sense because, in Gj, if a vertex n G 14 is in 
a component Gj then all vertices of I4 are in Gj. Consider the graph Gj together with all 
vertices and edges removed in the previous j iterations, i.e., the current Gq. Then, there exists 
some i, such that all cycles in Go are contained in the subgraph of G induced by the vertices 
of Gj. We claim that in Go, any edge that can be added without violating the upper bounds 
given by D, must have both its endpoints in some C(Gj). 

To prove the latter claim, we show that all other possible cases fail to happen. Assume 
that there exists an edge uv that can be added to Go without violating any constraint, such 
that u G V{Gi),v G C(Gj/), where i 7^ 4 . But then, uv would have been found and added 
in Case 1 of the current iteration. In particular, if I 4 G Mj and Vy G Ai', then we must 
have dxy = 0 or else Gj and Gj/ would become connected in Case 1 of the current iteration. 
Now, assume that there exists an edge uv that can be added to Go without violating any 
constraint, such that either u G V{Ci),v E Vy E Vi' with i' < j, or u E V^ E Vi,v E Vy E Vi> 
with i' < i < j. This means that a few iterations back, when G was Gj/, uv was available also. 
But Vy intersects at least two connected components of Gj/ and therefore, 3 n' G 14 such that 
u and v' he in different connected components in Gj/. Since adding uv is legal, so is adding 
uv'. But then, uv' would have been added in Case 1 of that iteration. This proves the claim. 

By the above claim, for any 14,14 fhat are not contained in V, we must have as many 
edges as possible between I4 and I4, that is, if G[X,Y] denotes the bipartite graph induced 
by X and Y , ^ 

\E{G0[V^,Vy])\=d^y. 
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Also, for any 14 and Ai we must already have as many edges as possible between 14 and all 
the 14 ’s in Ai in Gq- That is, 

|S(G„[y,U,^i,^,^yJ)|=max{|i4|, ^ 1 ,}. 

y:VyG^i 

Now notice that if we identify all Ai with one single vertex to get H from Go, H remains 
disconnected, but contains no cycles. That is, \E{H)\ < \V{H) \ — 1. But, 

|V'(//)| = |.4| + ^ \Vl and 

i-AiT 

1 _ _ ^ ^ ^ ^ 

x-y^^^T y-Vy^T x-.%^T V-VyGAi 

That is, if we use E and A to construct the weighted graph Q = (V,T,tc) as in the necessary 
and sufficient condition, then 


5^m(e)<|.4|+ Y. Il/.l-l, 

i.e., the condition fails to hold and {E,A) is indeed a certihcate showing that no connected 
graph in (V, d, D) exists. 

From the above, within no more than 2k iterations. Algorithm Valid Tree Construction 
either reduces the number of components of Go, or terminates, giving a certificate {E,A). 
Therefore, in at most 2fc|V| iterations the algorithm terminates, and if a valid tree T of 
(V, D) exists. Go is going to be such a tree. □ 


4 A step towards sampling 

Towards uniform sampling from hi = (V, d, Zi)),we can dehne a natural Markov chain on 
n. Let G G n and u, u' be two vertices that belong to V for some i. Also let n, v' be 
vertices such that uv,u'v' G E{G), but uv',u'v ^ E{G). Then, the graph G' with E{G') = 
E{G) U {mu', u'v} \ {uv, u'v'} is still in hi. We call such an operation a legal switch and denote 
it as [uv^u'v'\uv',u'v\ (Figure]^. Two legal switches are distinct, if they produce different 
graphs. 

We dehne the Markov chain Ai as follows. Given a state/graph Xt = G E fl, hrst calculate 
the number of distinct legal switches, i{G). Notice that this can be done in O(n^) time. 

• W.p. 1/2, let Xt+i = G. 

• W.p. 1/2, choose one of the i{G) distinct legal switches. Perform the switch and let 

G' be the resulting graph. With probability let W+i = G', otherwise let 

At+i = G. 
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Clearly, J\4 is aperiodic. Notice that if vr is the uniform distribution on hi, for all x, ?/ G 
Q,7i{x)P{x,y) = 7i{y)P{y,x). That is, if Ai is irreducible, then its unique stationary distri¬ 
bution is the uniform distribution tt. 

Below we show that A4 is irreducible. Usually this is a trivial step, but here it is a bit 
more involved. We should note here that [32] independently proposed an alternative proof for 
the irreducibility of Ai. This proof however was flawed, as noted in [101, where an alternative 
proof is given. Our approach below is simpler, although not completely straightforward. Let 
Go and Gi be two arbitrary instances in O. We want to show that there exists a sequence 
of legal switches that applied to Gq gives Gi. First, we need to introduce some notation and 
terminology. In what follows X = G © Gi is the symmetric difference of the current graph 
G and Gi; initially, G = Gq. We will refer to the edges of E{G) \ E{Gi) as straight, and the 
edges of E{Gi) \ E{G) as squiggly. The graph X is the union of straight and squiggly edges. 
Also, we will refer to the edges of E{G) nU(Gi) as dashed, and to the rest of the edges (edges 
in neither G nor Gi) as dotted. Obviously, the depiction of the edges will reflect their names. 

If there exist vertices x,u,v of X such that xu is straight, xv is squiggly and u,v belong 
to the same Vi, then we call x a pairing node (Figure]^. 



Figure 1: A legal switch. Figure 2: A pairing node. 

Finally, notice that Vu G X the number of straight and the number of squiggly edges 
adjacent to v are equal, i.e., they are degx{v)/2. 

Lemma 7. The Markov chain Ai defined above is irreducible. 

Proof. By induction on |£'(X)|. If |U(X)| = 4, then X has to contain an alternating cycle 
of length 4, with straight and squiggly edges, with two non adjacent vertices in the same U- 
Thus we can go from G to Gi in one legal switch. 

Now assume that |U(X)| = 2/c > 4. 

Case 1: X contains a pairing node x. 

Subcase la: There exists a straight neighbor of u, say w, that is also a squiggly neighbor of v 
(Figure]^). Then, by switching [vx,uw\vw,ux] in G, we reduce \E{X)\ to 2k — 4. Thus, by 
induction, we can go from G to Gi by performing a sequence of legal switches. 

Subcase lb: There exists a straight neighbor of u, say w, that is also a dotted neighbor of v, 
or there exists a dashed neighbor of u, w, that is also a squiggly neighbor of v (Figure]^). 
Then, by switching [vx,uw\vw,ux] in G, we reduce \E{X)\ to 2k —2. Thus, by induction, we 
can go from G to Gi by performing a sequence of legal switches. 

Subcase Ic: There exists a straight neighbor of u, say w, that is also a dashed neighbor of v, or 
there exists a dotted neighbor of u, w, that is also a squiggly neighbor of v (Figure [^). Notice 
that by switching [ux, vw\uw, vx] in Gi, we get a graph Gj, such that \E{G © Gj^)| = 2k — 2. 
Thus, by induction, we can go from G to G[ by performing a sequence of legal switches. Then, 
by switching [ux,vw\uw,vx] in Gj, we get Gi. 
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Figure 3: Case 1. 


Now, we claim that always one of the above subcases holds. Assume not. That is, (i) 
all straight neighbors of u {y excluded) are also straight neighbors of v and (ii) all squiggly 
neighbors of v {u excluded) are also squiggly neighbors of u. Since, a; is a squiggly neighbor 
of u and a straight neighbor of n, (i) implies that degx{u)/2 < degx{v)/2 (Figure]^) and 
(ii) implies that degx{u)/2 > degx{v)/2 (Figure]^), producing a contradiction. 




Figure 4: Subcases la, lb, and Ic exhaust all possibilities for Case 1. 

Case 2: X contains no pairing node. 

Then, there exist two edges, one squiggly xu and one straight yv such that u, v, x, y are all 
distinct and x,y G Vi, u,v E Vj for some i,j (Figure]^. Notice that xv can only be either 



Figure 5: If no pairing node exists. 

dashed or dotted. Assume it is dotted. {The case where xv is dashed is symmetric.) Then, v 
is a neighbor of y in G, but not of x. Since x and y have the same degree in G, there exists 
some w such that xw is in E{G), while yw is not. 

Subcase 2a: The edge xw is straight and yw is dotted. Then, by switching [xw,yv\xv,yw] in 
G, we keep \E{X)\ = 2k and we create the pairing node x (Figure [^). By Case 1 above and 
induction, we can go from G to Gi by performing a sequence of legal switches. 
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Subcase 2b: The edge xw is dashed and yw is squiggly. Like above, by switching [xw, yv\xv, yw] 
in G, we keep \E(X)\ = 2k and we create the pairing node x (Figure [^). By Case 1 above 
and induction, we can go from G to Gi by performing a sequence of legal switches. 

Subcase 2c: The edge xw is dashed and yw is dotted. Then, by switching [xw,yv\xv,yw] in 
G, we have \E{X)\ = 2k + 2 and we create two pairing nodes x, w (Figure |^). By examining 
the pairing node x with u, v like we did in subcases la, lb, Ic above, we either reduce \E{X)\ 
to — 2 and proceed with induction, or we reduce \E{X)\ to 2k, while X has the pairing 
node w, in which case, by Case 1 above and induction, we can go from G to Gi by performing 
a sequence of legal switches. 




Notice that in Case 2, it is not possible to have xw to be straight and yw to be squiggly. 

□ 


5 Further Directions 

For the original degree sequence graphic realization problem, when a graphic realization of 
a degree sequence di > d 2 > ■ ■ ■ > dn exists, several interesting generalizations can be 
solved efficiently. For example, if there are costs on edges, then we can find a minimum cost 
realization and we can generate a realization, or a connected realization uniformly at random 
(under mild restrictions on the degree sequence [22l [T9|, lU [T7]). 

The above problems have efficient algorithms because the original problem has a reduction 
to perfect matchings [221 ESI El]- In particular, for di > d 2 >...> dn as above, define the 
graph Q = (V,T) as follows. The vertex set is V = where Vi = ^j£[n],ji^i{vij} 

{uij/}, Wl < i < n. The vertices Vij denote a potential edge between i and j in the graphic 
realization. The vertices will enforce the required degrees di. Now the edges are £ = 
l-ii<i<j<n{{vij,Vji}} {{vij,Uij'}}. It is straightforward to verify that the 

degree sequence is realizable if and only if Q has a perfect matching. 

Is there a reduction from the joint-degree matrix realization problem to some version of 
matching, flow, or a similar better understood combinatorial problem? The relatively smooth 
decision and construction algorithm outlined in Section for realizability suggests that such a 
reduction might exist. One should probably try to reverse engineer the construction algorithm; 
however, the main difficulty is outlined in Remark 1 at the end of Section namely, the 
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alternating seqnences of edges involved in the algorithm are not pnre augmenting paths. The 
existence of a such a reduction may help solve the problems listed in the next paragraph. 

The most interesting open question is undoubtedly whether one can efficiently sample from 
(V, d, D). It would be interesting to prove that the Markov chain we suggest in Section or 
any other Markov chain for that matter, is rapid mixing. Although some progress has been 
made in this direction in [15], the problem remains largely open. Similar questions can be 
asked for a weighted version of the problem. Let (V, d, D) be an instance of the joint-degree 
matrix realization problem. If there is a cost associated with every potential edge, can we 
construct a realization of minimum cost? More importantly, can we generate uniformly at 
random such a realization? 

It is also natural to dehne the following generalization of the graphic realization problem, 
where conditions on the number of edges involve arbitrary subsets of vertices. In particular, 
for V = [n] and positive integers di > d 2 > ... > d„, let S = {Si, S 2 , ■ ■ ■, Sk} be an arbitrary 
partition of V, and let D = {dij), l<i<j<k,hesikxk matrix, where dij is the number 
of edges between Si and Sj. As is, this is the same as the partition adjacency matrix problem 
introduced in [9], and the skeleton graph problem introduced in [13], where some special cases 
are studied. Is there a polynomial decision/construction algorithm for the above problem in 
general? Notice that this is not a direct generalization of the joint-degree matrix problem, as 
was noted in [TH], unless we modify the above dehnition to specify the degree subsequences 
in each Si. 


References 

[ 1 ] W. Aiello, F. Chung, and L. Lu. Random evolution in massive graphs. FOCS, 2000. 

[2] W. Aiello, F. Chung, and L. Lu. A random graph model for massive graphs. STOC, 

2000 . 

[3] D. Alderson, J. Doyle, R. Govindan, and W. Willinger. Toward an optimization-driven 
framework for designing and generating realistic internet topologies. ACM Sigcomm CCR, 
2003. 

[4] A. Barabasi. Linked: The New Science of Networks. NY Perseus, 2002. 

[5] C. Berge. Graphs and Hypergraphs. North Holland Publishing Company, 1973. 

[6] K. Calvert, M. Doar, and E.W. Zegura. Modeling internet topology. IEEE Communica¬ 
tions Magazine, 1997. 

[7] F. Chung and L. Lu. Complex Graphs and Networks. AMS, 2004. 

[8] C. Cooper, M. Dyer, and C. Greenhill. Sampling regular graphs and a p2p network. 
SODA, 2005. 

[9] E. Czabarka. On realizations of a Partition Adjacency Matrix. Unpublished Manuscript, 
2014. 

[10] E. Czabarka, A. Dutle, P. L. Erdos, and 1. Miklos. On realizations of a joint degree 
matrix. Discrete Applied Mathematics, 181:283 - 288, 2015. 


16 


[11] S. N. Dorogovetisev and J. F. F. Mendes. Evolution of Networks: from Biological Nets 
to the Internet and WWW. Oxford University Press, 2003. 

[12] R. Durrett. Random Graph Dynamics. Cambridge U. Press, 2006. 

[13] P. L. Erdos, S. G. Hartke, L. van lersel, and L Miklos. Graph realizations constrained by 
skeleton graphs. ArXiv e-prints: 1508.00542, 2015. 

[14] P. Erdos and T. Gallai. Graphs with prescribed degrees of vertices. Math. Lapok, 11, 
1960. 

[15] P. L. Erdos, I. Miklos, and Z. Toroczkai. A decomposition based proof for fast mixing 
of a markov chain over balanced realizations of a joint degree matrix. SIAM Journal on 
Discrete Mathematics, 29(1):481-499, 2015. 

[16] M. Falontsos, P. Falontsos, and G. Falontsos. On power-law relationships of the internet 
topology. Sigcomm, 1999. 

[17] T. Feder, A. Gnetz, M. Mihail, and A. Saberi. A local exchange markov chain on graphs 
with given degrees and applications in p2p networks. FOGS, 2006. 

[18] M. Gjoka, B. Tillman, and A. Markoponlon. Gonstruction of simple graphs with a target 
joint degree matrix and beyond. In 2015 IEEE Gonference on Gomputer Gommunications 
(INFOGOM), pages 1553-1561, 2015. 

[19] G. Gkantsidis, M. Mihail, and E. Zegura. The markov chain simulation method for 
generating connected power law random graphs. ALENEX (SODA Workshop), 2003. 

[20] S.L. Hakimi. On the realizability of a set of integers as degrees of the vertices of a graph. 
SIAM, J. Applied Math., 10, 1960. 

[21] V. Havel. A remark on the existence of hnite graphs. Kaposis Pest Mat, 80, 1955. 

[22] M. Jerrum and A. Sinclair. Fast uniform generation of regural graphs. TGS, 73, 1990. 

[23] G. Jin, Q. Ghen, and S. Jamin. Inet: Internet topology generator, 

http://irl.eecs.umich.edu/jamin, University of Michigan technical Report, GSE-TR-433- 

00 , 2000 . 

[24] L. Li, D. Alderson, W. Willinger, and J. Doyle. A hrst-principles approach to under¬ 
standing the internet’s router level topology. Sigcomm, 2004. 

[25] P. Mahadevan, D. Krioukov, K. Fall, and A. Vahdat. Systematic topology analysis and 
generation using degree correlations. Sigcomm, 2006. 

[26] P. Mahadevan, D. Krioukov, M. Fomenkov, B. Huffaker, X. Dimitropoulos, kc daffy, and 
A. Vahdat. The internet as-topology: Three data sources and a dehnitive metric. AGM 
Sigcomm GGR, 2006. 

[27] A. Medina, I. Matta, and J. Byers. On the origin of power laws in internet topologies. 
AGM Sigcomm GGR, 2000. 


17 


[28] M. Newman. Assortative mixing in networks. Phys. Rev. Lett. 89, 2002. 

[29] M. Newman. Mixing patterns in networks. Phys. Rev. Lett. E 67, 2003. 

[30] M. Newman, Barabasi A.L., and D. Watts. The Strueture and Dynamies of Networks. 
Princeton U. Press, 2006. 

[31] A. Sinclair. Algorithms for Random Generation and Counting: A Markov Chain Ap¬ 
proach. Springer-Verlag, 1997. 

[32] I. Stanton and A. Pinar. Constructing and sampling graphs with a prescribed joint degree 
distribution. J. Exp. Algorithmics, 17:3.5:3.1-3.5:3.25, September 2012. 

[33] L. Valiant. The complexity of enumeration and reliability problems. SIAM, J. Computing, 
8 , 1979. 

[34] E. W. Zegura, Ken Calvert, and S. Bhattacharjee. How to model an internetwork. In 
Infocom, 1996. 


18 



